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PREFACE 


This report details the control analysis phase of the project 
"Modeling and Control of Transonic Cryogenic Wind Tunnels," sponsored by 
NASA/Langley Research Center (LaRC) under research grant NSG-1503. This 
work was performed during the project period ending March 1980. The 
contents of this docviment complement the modeling phase activity which 
has. been reported as "Synthesis of a Control Model for a Liquid Nitrogen 
Cooled Closed Circuit Cryogenic Nitrogen Wind Tunnel and its Validation" 
(ref. 4). This document reports the details of control law design, proof 
of its adequacy, microprocessor compatible software design, and electronic 
hardware realization and its successful performance on the 0.3-m 
Transonic Cryogenic Tunnel at NASA/LaRC. 
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AUTOMATIC CONTROL OF A LIQUID NITROGEN COOLED, 
CLOSED-CIRCUIT, CRYOGENIC PRESSURE TUNNEL 

By ^ 

S, Balakrishna^ 

INTRODUCTION 

Need for realizing full flight Reynolds number airflow past scaled 
models in transonic wind tunnels has been keenly felt for some time now, 
particularly to aid in understanding shock-boundary layer interaction. 

In the 1970 *s., there was considerable review of various methods for 
obtaining high Reynolds number airflow in relatively small wind tunnels 
without any fan power or dynamic pressure penalties. The most promising 
among many possible techniques proved to be the concept of operating 
the test gas of a wind tunnel at cryogenic temperatures (ref. 1). This 
was elegantly demonstrated by the proof-of-concept operation of the 
0,3-m Transonic Cryogenic Tunnel (TCT) facility at NASA/Langley Research 
Center (LaRC) in 1974 (ref. 2). Nitrogen was used as the test medium 
and was proven to behave perfectly under isentropic flow conditions 
(ref, 3). Further, the successful operation of the 0.3-m TCT was aided 
by contemporary advances in dryogenic materials and instrtimentation 
techniques. 

The initial operating experience on the 0J.3-m TCT facility indicated 
that the dryogenic tunnel process was basically well behaved from a 
control point of view, but exhibited continuous creep of the tunnel 
variables. This was caused by unmatched mass-energy control inputs which 
integrated to large errors in time when controlled manually. Further, 
difficulties were experienced in manual control for large changes in the 
tunnel conditions because of considerable cross-coupling between the 
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tunnel flow variables and the control inputs, making the tunnel runs 
rather uneconomical. Hence, the need for design of a fast, automatic, 
closed-loop controller for regulating the tunnel variables was keenly- 
felt . 


For these reasons, studies have been undertaken to formulate and solve 
the problem of closed-loop control of a fan-driven, liquid nitrogen 
cooled, closed-circuit, cryogenic pressure tunnel, with specific reference 
to the 0.3-m TCT facility at NASA/LaRC. This problem has been addressed 
in two phases: first, by the synthesis of a control compatible mathematical 

model of a cryogenic pressure tunnel starting from the basic physical 
laws and its experimental validation. The second phase was the design of 
a closed-loop, feedback controller to regulate the tunnel.^ s total tempera- 
ture and total pressure, the controller hardware realization, and its 
performance on the 0,3-m TCT facility. 

The effor' 3 of the synthesis of a control model of cryogenic tunnel 
process and its experimental validation have been successfully completed 
and reported in reference 4. This synthesis was based on the thermody- 
namic analysis of the tunnel resident gas and metal mass and their interaction 
with the various control inputs. The resultant lumped parameter, non- 
linear, multivariable model was used to develop a cryogenic tunnel simulator 
on a hybrid analog/ digital computer. The synthesized mathematical model 
was validated by reconciling the quasi-steady state and transient impulse 
response of the 0,3-m TCT facility and its real-time simulator responses 
under similar operating conditions. 

In this report, the complementary problem of the control system design, 
performance analysis for its accuracy and stability, microprocessor-based 
controller software development, the microprocessor hardware specifications, 
and finally the performance of the designed controllers on the 0.3-m TCT 
facility have been detailed. Firstly, the control model has been used to 
design two single- input single-output (SISO) controllers for the automatic 
control of the tunnel's total temperature and total pressure by assuming 
small perturbation local linearity. In view of the global nonlinear behavior 
of the tunnel model, the SISO control analysis has resulted in nonlinear ly 
gain scheduled proportional- integral-derivative (pID) control laws for 
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each of the temperature and pressure control loops. Strong cross-coupling 
of the fan power into the temperature loop has been countered by a fan 
power feed-forward feature. Further, in order to allow for directional 
nonlinearity, use of error magnitude dependent logic which switches the 
feed- forward on or off as necessary has been incorporated. Appropriate 
stability and accuracy tests have been performed to assure desired 
performance of the analytically designed control laws. 

The control problem, has also been addressed as a linear multivariable 
control problem, again by assuming small perturbation local linearity. A 
noninteractive controller so designed is shown to embrace the SISO control 
laws as a particular case. 

These control laws have been tested on the real-time, hybrid, cryogenic 
tunnel simulator for the adequacy of their performance. The SISO laws 
have demonstrated their ability to hold the tunnel total temperature to 
±0.25 K and total pressure to ±0.017 atm, with a typical step response 
settling time of about 20 sec for small amplitude step commands. 

Further, the SISO control laws have been translated to a microprocessor- 
compatible, functional flow chart which is presented in Appendix A. The 
microprocessor hardware specifications have been generated by looking at 
the electronic system design as a direct digital controller. The operator 
control interface design has been detailed. Finally, the successful 
operation of the microprocessor-based tunnel closed-loop controller on the 
0.3-m TCT facility, throughout the tunnel operational envelope, has been 
demonstrated. The ability of the controller to hold ±0.25 K in temperature 
and ±0.017 atm in pressure has been shown. 

The analysis and design of the tunnel test section Mach number automatic 
control loop is presented in Appendix B. The existing fan speed control 
system has some singularities in fan speed wherein the fan cannot be operated. 
Though a design of a Mach loop controller has been presented, it cannot 
be implemented till appropriate improvements are made to the fan speed 
control system. 
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NOMENCLATURE 


a 

XX 


A 

b 

b 


XX 


C 

P 

G 



GN2 

h 

H 

J 


elements of matrix Z 
area (m^ or m^/m^) 
tunnel circuit loss factor 
elements of matrix 

specific heats of nitrogen gas (KJ/kg-K) 

feed-forward path model 

gaseous nitrogen 

enthalpy (KJ/kg) 

feedback path model 

Joules ki lojQiile^X-. 


K,k 

LN2 

« 

m 

M 

N 

P 

Q 

5 
T 
t 
U 

V 

w 

Z ,R,C 
a 

6 

Y 
A 


gain constants (K also stands for Kelvin) 

liquid nitrogen 

mass flow rate (kg/sec) 

test section Mach number 

fan speed (rpm) 

pressure (atm) 

heat flow rate (KJ/sec) 

Laplace operator 

temperature (K) 

time or time constant (sec) 

internal energy (KJ) 

volume (m^) 

mass (kg) 

system matrices 

cooling capacity of gaseous nitrogen (KJ/kg) 
cooling capacity of liquid nitrogen (lU/kg) 
ratio of specific heats 
increment 

transport lag (sec) 



0 

thermal mass (KJ/R) 


i 

integer index 


• 

dot over a quantity refers to of that quantity 


r 

pressure ratio 


n 

fan efficiency 


a + j oj 

complex plane point; j 


y 

micron 


Controller 


T 

tunnel total temperature (K) 


ST 

set point, temperature loop (K) 


ET 

error, temperature loop (K) 


pT 

proportional gain, temperature loop 


DT 

derivative gain, temperature loop 


IT 

integral gain, temperature loop 


KT 

temperature loop gain 


QDF 

fan power (KJ/sec) 


KL 

full open liquid nitrogen flow (kg/sec) 


ALQ 

liquid valve area 


UT 

temperature loop update rate 


PL 

liquid nitrogen pressure (atm) 


P 

tunnel total pressure (atm) 


PS' 

tunnel test section static pressure (atm) 


SP 

set point, pressure loop (atm) 


EP 

error, pressure loop (atm) 


PP 

proportional gain, pressure loop 


DP 

derivative gain, pressure loop 


IP 

integral gain, pressure loop 


KP 

pressure loop gain. 


MAX 

maximum allowed liquid nitrogen valve area (percent) 
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KB 

bias gain 

FF 

feed-forward binary logic, on- 

LLT 

liquid limit binary logic, on- 

AG 

gas bleed valve area 

UP 

pressure loop update rate 

UA 

analog valve update rate 

KLl 

liquid valve gain, KL = KLl ^ 

KG 

gaseous valve gain 

KF 

fan gain 

TGR 

temperature gradient (K/sec) 

SISO 

single input, single output 

SM 

set point, Mach loop 

EM 

error, Mach loop 

KM 

loop gain, Mach loop 

IM 

integral gain, Mach loop 

pM 

proportional gain, Mach loop 

NS 

set point, fan speed loop 


Subscripts and Superscripts 
F,f fan 

g gas 

L liquid 

ra metal or motor constant 

s static 

1,2 analog bleed valves, or sensors 

* test section 

(n) nth sample where n is an integer = 0,1,.,. etc. 

set setpoint 

past previous valve 
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CRYOGENIC TUNNEL PROCESS 


A closed-circuit wind tunnel is a device intended for testing scaled 
models in fluid flow in order to obtain their aerodynamic characterictics* 

The test section fluid flow over the model is associated with certain flow 
parameters to which the aerodynamic characteristics of the model are related. 
For meaningful scaling of model data to its full scale, it is necessary 
to maintain the flow similarity between the wind-tunnel test section flow 
and the full-scale flow. Among many flow parameters necessary for flow 
similarity, Mach number and Reynolds number are the most important flow 
parameters for stationary models. For the first time, full-scale flow 
Reynolds number can be obtained at transonic speeds in small. wind tunnels 
because of the advent of the cryogenic flow wind tunnel concept (ref. 1). 

The operation of a closed-circuit, fan-driven, transonic, cryogenic 
nitrogen tunnel and its ability to allow independent control of flow 
Mach number, Reynolds number, and flow dynamic pressure have been proven 
by the operation of the 0.3-m TCT facility at NASA/LaRC shown in figure 1 
(refs. 1 and 2). This has been possible because of the ability of cryogenic 
tunnels to allow independent control of the tunnel gas temperature, the 
tunnel total pressure, and the test section mass flow achieved by three 
independent control inputs: viz liquid nitrogen flow into the tunnel, 

gaseous bleed out of the tunnel, and the fan speed. 

A closed-circuit, cryogenic wind tunnel basically consists of a 
thermally autonomous pressure vessel designed as an aerodynamically contoured 
endless duct in which the fan operation results in a range of flow 
velocities at the test section. The thermal autonomy of a cryogenic tunnel 
is realized by thermal insulation either external to the tunnel metal 
walls or internal to the metal walls. The operation of the fan in the 
tunnel results in motion of the test gas with the accompanying compression 
heating and wall friction heating. In order to maintain the gas temperature 
in the presence of this heating, the liquified form of the test gas, liquid 
nitrogen (LN 2 ) , is sprayed into the tunnel. The liquid nitrogen mass flow 
control allows either a cool down or regulation of the* tunnel gas tempera- 
ture. However, in this cooling process, the mass of the tunnel resident gas 
increases. In order to maintain the tunnel pressure, it is necessary to 
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bleed the warmer tunnel gas out of the tunnel, and this, controlled gas 
bleed mass flow allows tunnel pressure control. This complex interaction 
of the fan- induced heating of the tunnel gas, the liquid nitrogen inflow 
induced enthalpy decrease and mass increase, and the gaseous nitrogen 
bleed induced mass and enthalpy decrease, constitutes the tunnel control 
problem. 

This thermodynamic mass-energy interaction has been analyzed in great 
detail in reference 4, by graphical and analytical techniques. Such an 
analysis has resulted in a lumped multivariable model of a cryogenic 
tunnel. The role of the dominant energy storage term, the metal wall 
stored enthalpy, has been shown to be significant. The following equa- 
tions, drawn from reference 4, present the basic behavior of a closed- 
circuit pressure tunnel. 

Any closed-circuit wind tunnel with a finite pressure volume V, 
with nitrogen as the test gas, has a resident gas mass of 

Wg = 338.9 Y ^ 


when the tunnel fan is run up to a speed of N rpm, the fan develops a 
temperature and pressure rise across itself by consuming a finite amount 
of power. Considering the fan compression to be isentropic, we have 


temperature rise AT 



Jcelvin 


where 


( 2 ) 


P 

fan outlet 
^fan inlet 
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This fan pressure ratio for a closed-circuit tunnel is a function of 
tunnel circuit losses and can be expressed as 

pressure ratio r = 1 + bM^ (3) 

where for the 0,3-m TCT facility, the experimentally determined steady- 
state^loss factor is 


b = 0.197 (l -^) 

This fan pressure ratio term is a function of Reynolds number and has been 
expressed empirically in the above equation. 

The fan consumes a finite power Q in creating the tunnel flow, 

r 

which can be expressed as 

Q = KJ/sec (4) 

(1 + 0.2m2)3 


where 




6965 * 

A 

n 


C b 
P 



* 

A = test section throat area 


Further, the fan speed and the test section mass flow are related by an 
expression 

N = k Vt M (5) 

m ” 

where, for the 0.3-m TCT,, the experimentally derived constant is, 


-0.035 

k = 597(1 - 0.3M)P 
m 
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In order to maintain the tunnel gas temperature, a liquid nitrogen 
source of pressure is used to generate a controlled liquid nitrogen 

mass flow of m^ into the tunnel. This is obtained by varying the liquid 
nitrogen flow control valve area A^. 




( 6 ) 


This incoming liquid mass flow is associated with a specific cooling 
capacity 6, which is the heat absorbed by unit mass of liquid nitrogen 
when taken from the source and evaporated to tunnel temperature at fixed 
tunnel *^)ressure (ref. 5). Then for equilibrium conditions we have 


“l 6 ‘ Op • '’t 


(7) 


where 


6 = ^121 + CpT^ KJ/kg 


( 8 ) 


and 


= Heat absorbed by metal wall 

In order to maintain the tunnel pressure, the warmer tunnel gas is 

bled out at a mass flow rate of m by a throttle valve of controlled area 

g 

A , out to atmosphere. This mass bleed has a finite enthalpy loss and 
g 

has a cooling capacity of a. When the tunnel is in equilibrium we have 



(9) 


where, for the.0,3-m TCT, the bleed valve has choked flow K = 21.8 for 

g 

P > 1.5 atm, and unchoked flow 


= 21.8 2 - (^) for 1 < P < 1.5 atm 
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The cooling capacity of the gas bleed is 


■ (% - '=v) " 


== 0.3 T - 0.08 P 


lU/kg 
ICJ/kg j 


( 10 ) 


For the 0.3-m TCT, liquid nitrogen and gas bleed into and out of the 
tunnel are managed as shown in the schematic figure 2'. 

In view of the endless nature of the tunnel duct, the total tunnel 

resident gas mass W crosses a given section of tunnel at the rate of 

8 . 

the tunnel test section mass flow m. Thus, under a one- dimensional flow 
situation, a particle moves around the circuit in a finite time called 
circuit time under steady-state conditions. 


0.0486 V 


W 

t = -^ = * 

m A m/t* 


1 + 250 — 

1 + 0.2 m2 

t2. 



sec 


( 11 ) 


The estimated tunnel circuit time for the 0.3-m TCT facility is shown in 
figure 3, and varies as a function of the tunnel conditions. The fastest 
circuit time of 0.55 sec occurs at M = 1, P = 1 atm, and T =300 K. The 
slowest circuit time of 3 sec occurs atM=0.2, P=5 atm, and T = 100 K. 

The tunnel resident metal, which is exposed to the gas flow, has a 
finite amount of enthalpy which is released into the gas whenever the gas 
temperature differs from the metal temperature. This heat transfer is 
not uniform and occurs throughout the internal surface of the metal in the 
tunnel. The mode of heat transfer is the dominantly forced convection 
type. This heat flow mechanism has been detailed in reference 4, and can 
be expressed as a simple first order model: 

W C TS 

\ ■ (TtVS) Wsac (12) 

m 

where t = metal time constant. 
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For the 0.3-m TCT facility, this metal-gas. heat transfer time constant 

varies as a function of tunnel P, T, and M. A plot of t is shown 

m 

in figure 4. The largest time constant of 500 sec occurs when heat 
transfer is slowest at M = 0.2, P = 1 atm, and T = 300 K. The smallest 
time constant of 20 sec occurs when the heat transfer is best at M = 1, 

P = 5 atm, and T = 100 K. From this metal time constant, a gas time 
constant can be defined as 


W C 

t = — ■■ ^ ^ ■ t sec 
^g W C + W^C m ^ 
® g V t m 


(13) 


A plot of the estimated gas time constant is shown in figure 5 and corresponds 
to the 0.3-m TCT facility; The gas time constant is smallest at M = 1, 

T = 300 K, and P = 1 atm. The largest time constant occurs at M = 0.2, 

T = 100 K, and P = 5 atm. 

The cryogenic tunnel thermodynamic model, shown in figure 6, can now 

be used to derive the cryogenic . tunnel dynamic model. 


Energy 




(14) 


where 


= liquid nitrogen particle transit time to the big end 

== 0.8 t sec 
c 

= gas pulse particle transit time - 0.1 t^ sec 
= fan heat transit time ~ 0.6 t 

c 

= W C + W C 

g V t m 
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Mass 


dt 




( 15 ) 


where K^bMg corresponds to fan acceleration contribution to the pressure. 
Tunnel/Fari 


M 



-T S 

H e ^ 

ik (1 + t^S) 
in p- 


(16) 


where 


= acoustic time lag from fan to big end 

t = equivalent first order plenum time constant 

k = 597(1 - 0.3 m) 
m 

These basic equations provide the lumped equivalent of the tunnel 

dynamical behavior, by expressing tunnel variables T, M, and P as 

functions of the tunnel control inputs A- , N, and A . Obviously a 

^ a 

considerable amount of cross-coupling exists amongst the tunnel variables 
and the inputs. These equations have been put together to form a multi- 
variable matrix model of a cryogenic pressure tunnel, which is shown in 
figure 7, 

This model has seven nonzero elements and associated gain terms. 

Equation (12) describing the tunnel temperature dynamics has three gain 

terms. The steady-state temperature gain for LN 2 - input is- (0 + a )/e. 

A plot of this gain for the 0.3-m TCT is shown in figure 8. Typically, 

the steady-state sensitivity of the tunnel for LN 2 input is about -0,12 K/kg. 

To cool the tunnel from 300 K to 100 K, one needs about 1,666 kg of liquid 

nitrogen. However, under transient conditions, the tunnel temperature 

gain term is (B + ct)/W C . This transient gain for the 0.3-m TCT is 

S V 

shown in figure 9, and varies from -1 K/kg to -40 K/kg. The steady-state 
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temperature gain for gas bleed is -a/0 and is shown in figure 10. Since 
|a| < |3|, its contribution is very small. 

The tunnel Mach number steady-state gain is a function of K and 

CT ^ ^ 

1/^. It is not affected by liquid or gas inputs directly, but coupling 

does occur through 1 /^t". 

The cryogenic tunnel pressure gains are basically derived from 

equation (15) , and are functions of two opposing pressure gradient terms 

P • P • 

— A m and ^ • Together they provide the characteristic pressure 

8 

signature of a cryogenic tunnel for a liquid nitrogen input. This has been 
extensively detailed in reference 4. The tunnel pressure gain for gas 

p 

bleed control input is a simple function of — due to insignificant con- 

p . 8 

tributions from — T. Further, the tunnel pressure gain for fan speed 

^ P * dM 

input is contributed to by — T and bMp — . 
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CONTROL DEVICES 


In order to control a cryogenic pressure tunnel, it is necessary to 
have control devices which allow fast and linear control of the liquid 
nitrogen flow into the tunnel, gaseous nitrogen bleed mass flow out of the 
tunnel, and an adequate degree of, resolution and speed of response in the, 
fan speed control. Further, it is necessary to sense the tunnel variables — 
total pressure, total temperature, and test section Mach number, to a 
high degree of resolution, accuracy, and speed of response. The details 
of the control and sensing devices used for the 0.3-m TCT facility are 
now presented. 


Liquid Nitrogen Flow 

The operating envelope of the 0.3-m TCT covers a tunnel pressure of 
6 atm at 300 K and a test section Mach number of about 0,8. Under these 
conditions, the fan power consumption is very high and has been estimated 
from equation (4). To cancel this heat, a liquid nitrogen flow of 

. . / 

^ ^ kg/sec 

is necessary. The maximum flow is about 10 kg/sec. The liquid nitrogen 
mass flow is estimated from the equation 

Kp m3 

m = 

(1 + 0.2m2)3(121 + T) 


This has been plotted as a function of tunnel test section Mach number and 
is shown in figure 11, which represents the tunnel under equilibrium 
conditions. 

The liquid nitrogen injection system must hence be capable of delivering 
a mass flow of 10 kg/sec into a 6— atm pressure vessel. In the 0.3-m TCT 
facility, a positive displacement pump capable of 15 kg/sec flow at 10 atm 
has been used for the LN 2 system. This flow is derived from a 212,000-liter 
twin tank system, in which the tanks are double walled for vacuum insulation. 
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A schematic of this system is shown in figure 2. The pump pressure is 
maintained by a process control type relief valve driven pneiomatically . 

The drive signal is from a pneumatic pi controller, operating on the 
difference between the set point and the liquid pressure. This system 
has been tuned to 1 Hz bandwidth at critical damping. 

The liquid nitrogen flow is controlled by four digital valves. Each 
valve has a resolution of 1/1024 of full area. This is achieved by a 
set of 10 binary weighted orifices which are controlled by a set of 
solenoids. The flow orifice area starts at 1/1024 of full opening and 
progressively doubles in the subsequent 9 orifices. The liquid nitrogen 
mass flow, from the four digital valves which function in parallel, is 

where varies from 0 to 1 in steps of 1/1024. Each solenoid has a 

response time of about 20 msec, and they function in parallel. 

Gaseous Nitrogen Bleed 

Removal of the warmer tunnel gas from the settling chamber can be per- 
formed either on a passive mode or on an active mode. In the case of the 
0.3-m TCT, a passive bleed mode for removal of tunnel gas has been used. 
This is because the tunnel pressure is always in excess of one atmosphere. 
In order to maintain the tunnel pressure, the gas bleed valve must have 
the capability of removing a mass flow m which is generally greater 
than m^. Further, since the liquid valve has a resolution of 1/1024, 
the gas valve should possess somewhat similar resolution. 

As indicated in figure 2, the 0.3-m TCT facility has 3 pneumatically 
controlled valves operated through process control pneumatic controllers. 

In series with one of these valves, a 10— element, 8— bit digital valve has 
been used for purposes of automatic control. This valve has a resolution 
of 1 in 256. The sonic nozzles representing these binary weighted orifices 
are controlled in an on-off mode by solenoid piloted pneumatic actuators. 
This has been necessary because of the large nozzle areas needed to create 
gas mass flow at low pressure differentials. 
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Figure 12 shows the valve area necessary to maintain the tunnel at 
equilibrium. Obviously when the tunnel is running at low pressures, or 
under unsteady conditions, the valve area required exceeds 100 percent. 

It is under these conditions that the other two valves can be operated 
as bias valves, thereby enhancing the equivalent resolution of the digital 
valves to about 1 in 800. The valve area is estimated as 


A 

8 


TM^ 

(1 + 0.2m2)3 (121 + T) K 

8 




The response speed of the bigger elements of the digital valve is about 
150 msec. In a cryogenic tunnel, any valve with a resolution of about 
1 in 800 and a full area control response of about 0.5 sec can be 
considered adequate for pressure control. 


Fan Speed Control 

The closed-circuit wind tunnel test section Mach number is controlled 
by the mass flow, which in turn is controlled by the fan pressure ratio 
and the tunnel circuit loss characteristics. In the case of the 0.3-m 
TCT, the steady-state tunnel loss/fan characteristics have been experi- 
mentally derived and expressed as a fan speed-Mach relationship-^ 


H = M 

where 


k * 597 (1 - 0.3M) 
m 

A plot of the fan speed as a function of steady-state tunnel variables, total 
pressure, total temperature, and the test section Mach number is shown in 
figure 13. The fan has been designed to operate away from the surge line 
throughout the tunnel envelope. 

The tunnel fan speed control is achieved by controlling the electrical 
supply frequency to the two-pole induction motor driving the fan. In order 
to maintain motor magnetic flux density, the voltage/f requency ratio of this 
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variable frequency power generator is maintained constant. The speed is 
directly proportional to frequency as 120 x frequency /number of poles. 

The variable frequency supply is a rotating Kramer type of system. 

The tunnel variables: total pressure, total temperature and the 

test section Mach number are sensed by appropriate transducers of adequate 
accuracy and response. The 0.3-m TCT facility has the following tunnel 
variable sensing transducers. 

Tunnel Total Temperature 

The total temperature of the tunnel gas is sensed by a copper constantan 
junction thermocouple of 800-p wire diameter. This size is considered • 
a compromise between speed of response and strength of the sensor to 
withstand high dynamic pressure flows. This thermocouple is located at 
the center of the settling chamber and downstream of the flow-straightening 
screens. The output of the thermocouple is nonlinear. Hence the output of 
the thermocouple is referred to a reference junction maintained at ice point 
and is linearized to provide 0.1 percent linearity and infinite resolution. 
Secondly, the response time constant of the thermocouple is a function of 
the flow conditions over the temperature sensing junction head'. It is 
estimated that the response time constant varies from 0.5 to 3 sec over 
the tunnel operational envelope. 

Total Pressure 

The tunnel total pressure is sensed at the settling chamber, downstream 
of screens and upstream of the contraction using a total pressure probe. 

This probe is connected to a bonded strain gauge type pressure transducer 
through an instrument tubing. The transducer behaves linearly and, in 
the presence of the tube, has an estimated response time of 0.25 sec as 
a first order l^g for small perturbation, 

Mach Number 

The cryogenic tunnel test section Mach number is estimated on line 
using two pressure signals corresponding to the tunnel total pressure and 
the tunnel test section static pressure. The tunnel test section static 
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pressure probe is connected to a bonded strain gauge transducer through a 
tube and has a response time of 0.25 sec. Utilizing these pressures, 
the Mach nxomber is estimated as 



The Mach number accuracy is estimrted to be about 0.002 to 0.003, and 
the response about 0.25 sec or better. 

Figure 1 shows the locations of the control inputs, the location 
of the sensors, and figure 14 shows the actuator and sensor performance. 
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CONTROL ANALYSIS 


It is well known that classical and modern control theory provide 

the basis for the synthesis or analysis of controlling physical systems. 

In order to generate the appropriate control, strategies, it is necessary 

to describe the physical system in the form of a control-compatible, 

input-output type of mathematical model. The cryogenic tunnel has been 

so modeled in the previous sections, and it uses the three control 

inputs, liquid nitrogen valve area A^, fan speed N, and gaseous nitrogen 

bleed valve area A in order to control the cryogenic tunnel gas total 

S 

temperature, the test section Mach number, and the total pressure. A study 
of the cryogenic tunnel mathematical model of figure 7 reveals that the 
system is time invariant in parameters, multivariable, and nonlinear in 
its characteristics. Obviously the closed-loop control analysis of such 
a system is not a very straightforward one. 

In the present research effort, the nonlinearity in the tunnel model 
is ignored by considering the tunnel to behave linearly for small perturba- 
tion. The problem of closed-loop control analysis then becomes the design 
of a controller for a t ime- invariant , linear, multivariable, dynamical 
system. 


Two analyses are presented in this section, the first being a further 
simplification of the model by ignoring the cross-coupling and analyzing 
the control problem as three separate SISO control efforts. The tunnel 
total temperature T is assumed to be controlled by the liquid nitrogen flow 
valve area Aj^, the tunnel test section Mach number M to be controlled 
by the fan speed N, and the tunnel total pressure P to be controlled 
by the gaseous bleed flow valve area A . The second analysis is to seek 
a noninteractive multivariable controller. 


Single Input/Single Output Analysis 

If the cryogenic tunnel model of figure 7 were to be a truly diagonal 

matrix with the off-diagonal terms tending to zero, then each of the 

tunnel control variables T, M, and P could be uniquely controlled 

by the control inputs A^ , N, and A , respectively. However, since 

L g 
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the off-diagonal elements are not all zero, these can be considered as 
disturbance on SISO control loops. In the following analysis, the design 
of temperature and pressure control are considered on an SISO basis. 

Temperature control loop . - The cryogenic tunnel, total temperature, 
control loop schematic is shown in figure 15. In this scheme the tunnel 
temperature is controlled by: the liquid valve area A^. The cross-coupling 
from the other two inputs, fan speed N and. gas valve area A , are 

o 

shown as external disturbances. The tunnel total temperature T is 
sensed by a thermocouple providing linear, first order output of time 
constant t^. The process temperature and the set point are compared, and 
'the error in the temperature is determined. This error is gain scheduled, 
a design objective to be determined in this. analysis, and is operated on 
by a pID control law. The resultant signal is used to drive the liquid 
nitrogen valve, after conditioning the signal to account for valve sensitiv- 
ity which varies as a function of liquid nitrogen source pressure and 

tunnel pressure P at the injection point. 

Consider the tunnel temperature dynamics of equation (14) : 


dt 


1 + t S\ 
— ^ 

1 + t S / 
m / 





- K A 
g g -0 


-T 


e 


S 

g 



The tunnel temperature is controlled 
proportional to N) . Since | (6 + a) 
ignored. However, Q^,/© is dominant; 
bance is necessary. 

Let the temperature error (ET) 
loop gain 

liquid nitrogen valve gain 
pID controller terms 

transport delay 
sensor dynamics 
valve area 

forward path transfer function 


by A^, A^, and Q^, (which is 

I > |a|, the contribution of A is 

g 

hence a feed forward of this distur- 

= T - (ST) 

= (KT) 

= (KL) 

= (pT), (DT), (IT) 

. (1 - T^S> 

= 1/(1 + tiS) = H(S) 

= (A^) = ALQ 
= G(S) 
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The open-loop gain product of the temperature loop is 


G(S) H(S) = (KL)(A^) •• [(pT)(ET) + (IT) / (ET) + (DT) ^ (ET) 


(1 + t S) (1 - T,S) 

m ^ 

(1 + t S) ’ (1 + t,S) 
8 


» -(KK) 


j(DT) S2 + (pT) S + (IT)) jl + (t^ - T^) S - 


S2(l + tiS) (1 + t S) 
g 


(17) 


The system stability is determined by solving the polynomial G(S) H(S) + 
1=0 ana studying the root locations in the complex plane; 


0 = 1 + G(S) H(S) 

= S* [t,t^ - (KK)(DT) 

+ S3 [ti + tg - (KK)(DT)(t^ - T^) - (KK)(pT) t^xj 
+ $2 [l + (KK) j(DT) + (pT)(t^ - - (IT) 

+ S [(KK) (pT) + (KK) (IT) (t^ - Tj^)] 

+ (KK)(IT) 


where 


KK 


.-(KL)(A^)(S + a) 
0 


(18) 


The roots of the polynomial of equation (18), when (KK) ^ 0, are called the 
poles of the system. The poles are also the roots of the denominator poly- 
nomial of equation (17). The roots of the polynomial of equation (18), when 
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(KK) ^ are the zeros of the system and correspond to the roots of the 
numerator polynomial of equation (17) . As the system gain (KK) is taken 
from zero to a large value, the roots of the system start from the poles 
and end at the zeros. These root movements can be represented graphically 
on the complex plane, and these loci are called the root loci of the 
system. The location of the roots at any finite gain (KK) provide an 
insight into the stability and accuracy of the system. 

The poles of the. temperature. control loop of . the 0.3-m TCT cryogenic 
tunnel are 

S = 0;S = 0;S = --^;S = --^ 

- t tl 

•g ^ 

The zeros of the temperature control loop are 


c A. c Ji., c -(pT) ± V(pT)^ - 4 (DT)(TtT 

^"“t ~ 2(DT) 

m L 

The qualitative asymptotic root loci for the temperature control loop 
are shown in figure 16. The axes of the plot are the real and imaginary 
values of the root S. Typically, any point (S ~ a + ju)) on the plot 
corresponds to a time-dependent response of the type 


St 

e 


(a ± jo))t 
e 


(19) 


As long as a is negative with m at zero, the system response is over- 
damped convergence. A a which is negative with finite ±jtu corresponds 
to an oscillatory convergence. When -a = ±jo) the system is critically 
damped; when a = 0 with a finite jco the system is oscillatory and continues 
to oscillate at a steady value. When roots move to the right side with 
a positive, the system diverges. 

The temperature control system has four poles and four zeros (fig. 16). 
Hence the four root trajectories start at poles and end at zeros. The two 
integrator poles move from the origin to the negative real zero of metal 
and controller. These root values correspond to the time response of metal. 
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The other two roots corresponding to gas and sensor are on the negative 
real axis, but as the system gain is increased the two move towards each 
other and break away symmetrically into the complex plane. At low Mach 
numbers, these roots move to the right and cross over, suggesting oscilla- 
tory and unstable responses at high gains. When the tunnel flow Mach 
number is high, the roots move to the left, suggesting stable behavior. 

In both the cases, the roots move to the zeros and one of these moves 
around the Reimann sphere to arrive at either controller zero or transport 
zero. 

These root movements are shown in the asymptotic root locus plot 
of figure 16. Further, figure 17 shows one quadrant of the gas roots for 
various Mach numbers at T = 200 K and P = 1, 3, and 5 atm. The system 
tends to be unstable at high gains atM=0.2. AtM=0.4 and 0.9, the 

system is more stable. Figure 18 shows one quadrant motion of gas roots 
for varying temperatures. The following inferences can be drawn from 
figures 16, 17, and 18: 

(1) The temperature control loop tends to go unstable at high 
gains and low Mach numbers. 

(2) Tunnel pressure has relatively weaker effects on system 
stability and bandwidth. 

(3) Tunnel temperature control bandwidth improves with 
Mach number. 

(4) Tunnel temperature control is very stable at high Mach 
Numbers . 

A set of root loci was generated by varying T from 100 K to 300 K, 

P from 1 to 5 atm, and m from 0.2 to 1, to cover the bulk of the tunnel 
operational envelope for the 0.3-m TCT facility. The loop gain necessary 
to operate- the temperature control loop at critical damping was determined 
and a curve fit was made. The pID control gains were also scanned. Thus an 
acceptable nonlinear gain schedule was determined for (KK) as 
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since 



( 20 ) 


where 


(KT) = loop gain (j^) 

Hence, a relatively .large loop gain can be used when M - 1, P - 6 atm, 

A T = 100 K, and there is a small gain atM=0-2, P-1.5 atm, and 
T - 300 K. The range of gain change is about 1:20. 

The steady-state accuracy of the temperature loop is. very good because 
the system has two integrations as shown in the model of equation (17). 

The temperature loop also has a zero error for ramp command and hence 
possesses adequate accuracy of control. 

Pressure control loop . - The cryogenic tunnel, pressure control loop 
schematic is shown in figure 19. The tunnel total pressure is assumed to 
be controlled by gas valve area A . This control loop is disturbed by 
liquid nitrogen mass flow input and by the tunnel temperature gradient. 

The tunnel total pressure is sensed by a transducer which has a first order 
response with a time constant of t 2 seconds. The tunnel pressure is 
compared with the desired pressure set point and the error is determined. 
This error is gain scheduled, the gain schedule being the design objective, 
and is operated on by a pID control law. The resultant signal is used 
to drive the gas bleed valve after conditioning it for valve flow gain 
sensitivity. 

Consider the tunnel total pressure dynamics of equation (15): 


dp P dT ^ P • P • 

dt T dt W ™L ~ W “g 

g g 


+ 




bHpf 
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The tunnel pressure is controlled by A 

. dT . dM „ ® 

A-, — , and — . Hence 


dt’ 


directly and is disturbed by 


dt 


W 


K A + [Disturbances]] 


( 21 ) 


The tunnel pressure control behavior can now be analyzed. Let 


pressure loop error 
loop gain 

gas bleed valve gain 
pID controller terms 

transport delay 

pressure sensor 

forward path transfer 
function 

Then open-loop gain product is 


= (EP) = P ~ (SP) 

= (KP) 

= (KG) 

= (pP), (DP), (IP) 

-T S 

=e ^ -1 -tS;t - 0.1 t 

S o c 

= 1/(1 + taS) = H(s) 

= G(S) 


G(S) H(S) = — ' [(pPHEP) + (IP) / (EP) + (DP) ^ (EP)1 

W J? S2 dt J 

(1 - T S) (22) 

Is — 

1 H- t2S 

The stability of this loop in its automous condition can be analyzed by 
using root locus analysis. 

0=1+ G(S) H(S) = s3 [ta - (KK) (DP) T 1 

s * 

+ s^fl + (KK) (DP) - (D,P) (KK) T 1 

^ g-* 

+ .s[(KK) (pP) - (IP) (KK) Tg] 

+ (IP) (KK) (23) 
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where 


( kk ) . mxmil 

“gV? 


The poles of the pressure loop polynomial are 

S = 0; S = 0; S = - ~ 

^2 

The zeros of the pressure loop polynomial are 

, -(pP) ± J(pP)2 - 4(DP)(IP) 

q - . c - J 

^ T ^ ^ 2 (DP) 

g 

The asymptotic root trajectories for the pressure loop are shoxjn in 
figure 20. This root set has three trajectories as (KK) is varied from 
zero to a large value. Obviously, the tunnel pressure loop is stable 
but oscillatory because of dominant disturbances from liquid flow, tem- 
perature gradient, and fan effects. The two integration poles break 
away from the origin to the left as a complex pole pair and join together 
at high gains to reach the transport zero by moving through the Reimann 
sphere pole. The inferences that can be drawn from this analysis are 

(1) The pressure loop is basically stable, but always 
oscillatory with an under damped convergent response. 

(2) The tunnel flow Mach number does not affect the pressure 
loop directly, but does so through the temperature loop. 

(3) Dual integration in the loop provides high steady-state 
accuracy but always has an overshoot. 

A number of tunnel conditions were tried on the above analysis, and 
the pressure loop gains (KK) were determined at critical dumping for 
various tunnel conditions. A nonlinear gain schedule was determined for 
the pressure loop as 


(KK) oc 


1 
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since 


(KK) 


CKG) (AG)p ^ 


(AG) « 



p2 (KG) 


or 


(AG) 


(KP) 1 . 


where 


(KP) is the loop gain term. 


( 24 ) 


Multivariable Control Analysis 

Consider a linear time- invariant , multivariable system described by 
a square (n x n) transfer function matrix [Z(S)] which is nonsingular 
and hence invertible. The basic principle of deriving a noninteractive 
controller for such a system is by diagonalization. A controller matrix 
[C(s)] is derived such that the open-loop system matrix product [C(S)] 

* [Z(S)] becomes a diagonal transfer function matrix [R(S)], whose 
diagonal elements form the characteristic equations of the n single 
input-single output systems (ref. 6). 

Thus 

[C(S)][Z(S)]= [R(5)] 


or 


[C(S)1 = [2~1(S)][R(S)] 

Most of the frequency domain multivariable controller design problems exist 
because of the characteristics of the inverse [Z~KS)J and its properties. 
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occur 


T S 

For a typical system, prediction terms like pure lead or e 
in the [Z*^(S)] which are difficult to realize. These problems have 
led to analytical techniques like the concept of diagonal dominance 
(ref. 7) and analyses by inverse Nyquist array (ref. 8 ) in the design 
of multivariable controllers. 

In the present research effort, it has been assximed that [Z*“^S)] 
is physically realizable by ignoring the transport delay terms. 

If a matrix Z is a 3 x 3 matrix of type 


Z 


^11 ^12 ^13 
0 3-22 ^ 

^31 ^32 ^33 


then 


Z"1 = 


bi2 bi3 

0 b 22 ® 

t>31 ^32 b33 


1 

Det 


^ 22^33 ^ 13 ^ 32 -^ 12^33 -^ 22^13 


311 ^ 33 -^ 13^31 


-322331 3i2a3i-ai 1332 ana22 


where 


Det = a 22 [311333 - ai 3 a 3 i] 


For the 0.3-m TCT facility modeled by the equation of figure 7 , we have 
an: inverse matrix Z”^S) shown in figure 21 , in which contributions of 
the transport delay terms have been ignored. If the goal is to realize 
three independent loops for temperature, pressure, and Mach number described 
by [R(S)]: 


[R(s)] 


1 

S 

0 

0 


0 

1 

1 + t s 

p 

0 


0 

0 

1 

s 
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Then the controller [C(S)] is 
[C(S)] = [Z’Us)][R(S)] 

and the controller [C(s)] is described by figure 20, in which the 
differentiator S is absent in columns 1 and 3 and (1 + t S) is absent 

• P 

in column 2, 

It may be noted that, with this noninteractive controller, the 
three open-loop transfer functions become 



(1 + t S) 

m 

s2(l + t S)(l + t S) 
P g 



The element bn of the controller could be simplified to the SISO 
temperature control law of equation (20), and the element b 33 of the 
controller is the SISO control law of equation (24). Thus the noninterac 
tive controller embraces the SISO control laws as particular cases. 
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VALIDATION OF THE SISO CONTROL LAWS 


As detailed in reference 4, the mathematical model of the 0.3-ra TCT 
facility has been simulated on a hybrid computer. This real-time, 
interactive, cryogenic tunnel simulator is shown in figure 22. The hybrid 
scheme for simulation consists of obtaining the linear dynamical solutions 
on the analog computer which is updated by the nonlinear coefficients 
which are computed precisely on the digital processor. Figure 23 shows 
the basic scheme for hybrid simulation. The analog and the digital 
processors communicate with each other at 25 Hz. In figure 24, the 
temperature dynamics simulation is shown. Figure 25 shows the pressure, 
Mach number, and fan dynamics simulation. In figure 26 the valve control 
interface is shown. The SISO control laws are implemented on the hybrid 
tunnel simulator as shown in figure 27, and as detailed by equations 
(20) and (24). 

This cryogenic tunnel simulator has been validated for its steady- 
state, quasi-steady, and transient response by good reconcilation of 
the 0.3-m TCT responses and the simulator responses for cool down, warm- 
up, and open door impulse responses for all the^ control inputs (ref. 4). 

The tunnel temperature loop is closed as shown in figures 23 and 
27, and it includes the model for a thermocouple whose response time 
varies as a function of tunnel conditions. The tunnel pressure loop is 
closed as shown in figures 24 and 27 and includes the pressure sensor 
dynamics. Both the tunnel temperature and pressure loop worked 
satisfactorily. In figure 28, the tunnel total temperature and the 
total pressure are on closed-loop control with the fan speed controlled 
on a speed-control loop maintained at a constant speed. The initial 
tunnel conditions are P = 2.2 atm, T = 263 K, and M = 0.69. Firstly, a 
temperature set point command of +16 K has been imposed. The tunnel 
temperature has changed from 263 K to 279 K in about 20 sec. The gas 
temperature responds quickly. Since the tunnel metal is still near 
263 K, the gas temperature varies slowly from 278 to 279 K. During this 
period, the tunnel pressure loop opposes the mass variation and maintains 
the pressure at 2.2 atm. The pressure is disturbed in an oscillatory 
manner but recovers to the steady-state value in about 25 sec. The 
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tunnel Mach number tracks the tunnel temperature by the law After 

about 90 sec, the tunnel gas temperature set point is changed from 279 K 
to 263 K. The control system maintains the tunnel pressure, and the 
temperature settles at 263 K in about 25 sec. Since' the metal temperature 
is still near 263 K, no creep in gas temperature is evident. 

Secondly, the pressure loop: set point has been changed from 2.2 atm 
to 1.86 atm, while keeping fan speed constant and temperature set point 
at 263 K (fig. 28). The tunnel pressure smoothly drops to 1.86 atm and 
has a convergent oscillatory response as predicted. The pressure settles 
down in about 30 sec. The temperature loop fights the enthalpy change 
caused by gas valve and settles down when the pressure is steady. The 
tunnel Mach number is affected only by the temperature during the 
transition, and the change of pressure has very little effect on Mach 
number. Subsequently, the tunnel pressure set point has been changed 
back to 2.2 atm. It may be noted that the pressure response is slower 
in transition to a higher pressure. This is caused by the tunnel gas 
accumulation which occurs only at the rate of The temperature loop 

fights the coupling and settles down when the pressure settles, Mach 
number is also changed back to the starting value. The pressure coupling 
on Mach number is through the term P 0.^35^ which is very weak. 

Thirdly, the tunnel fan speed has been changed such that the Mach 
number varies from 0.69 to 0.59 while keeping the temperature loop set 
point at 263 K and pressure set point at 2. 2 . atm. It may be noted that 
the control loops maintain temperature to within ±0.25 K and pressure to 
±0.017 atm. The Mach number has been further changed from 0.79 to 0.59 
and from 0.59 to 0.69. In the presence of these gross power changes, 
the pressure and the temperature have been well regulated by the control- 
lers. 

In figure 29, a similar set of set point commands for temperature, 
pressure, and fan speed is shown. The tunnel temperature is at 90 K, 
pressure at 2.2 atm, and Mach number at 0.75. Again, the temperature 
and pressure control systems behave adequately. When the system 
temperature changes from 90 K to 82 K, pressure is changed from 2.2 atm 
to 1.86 atm, and Mach number from 0.75 to 0.85 through 0.65. 
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In summary, the control laws derived in equations (20) and (24) 
provide satisfactory closed-loop temperature and pressure control. These 
SISO control laws are adequate , to counter the strong cross-coupling between 
the tunnel process variables, and the control inputs. 
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TUNNEL CONTROLLERS 


Hitherto the detailed development of the control law for automatic 
control of the 0.3-m TCT facility was discussed. Further, the adequacy 
of these control laws to hold closed-loop control within desired speci- 
fications was proved by testing the control algorithms on the real-time 
cryogenic tunnel simulator. In this section, the mechanization of these 
control laws using electronic control hardware is discussed. 

Tunnel Closed-Loop Controls 

The detailed schematic diagrams for the tunnel temperature and 
pressure are shown in figures 30 and 31. 

Temperature loop. - In figure 30, the temperature control scheme used 
for the 0.3-m TCT facility is detailed. The control law used is 


(ET) = T-(ST) 


(ALQ) 


M . yp (KT ) I " ^ r ^ ^ 

(KL)Vt L 


+ (KB)(FF) 


(25) 


This control law is the one derived in equation (20). Additionally, provision 
has been made to have fan power bias on an on-off mode. The bias is the 
equivalent liquid nitrogen flow to cancel fan heat. Hence, 


(KB) 


This has been compensated for valve gain, as 


(kb) 

(121+T) (KL) 


(26) 


In a cryogenic tunnel, the fan power feed forward is ntcessary only while 
cooling down or while regulating the temperature. When the tunnel is warming 
up, feed forward would be a waste of liquid nitrogen. Hence the feed forward 
is switched on based on the following logic: 
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If 

(ET) > 0 

FF is on 


If 

(ET) < -5 

FF is off 

(27) 


This error--based switching of feed forward is shown in figure 32. This 
logic has a built-in 5 K hysteresis to prevent hunting in temperature. 

Secondly, when the tunnel is warm at 300 K and at fairly low pressure, 
switching the automatic temperature loop with a 100 K set point would result 
in large liquid nitrogen flow on a pID law. This may lead to liquid in the 
tunnel and large thermal shocks. It may be noted that to cool the 0.3-m 
TCT nearly 1,800 kg of LN2 is required and this cannot be done fast. Hence, 
a saturation term on the LN2 flow from the controller is proposed. 


Maximum LN2 allowed: 


Fan power ^ Steady gas temperature 

g (KL) cool down at a known gradient 


Max = (KB) (FF) (1-LLT) + 


255 PV (TGR) 
T(121+T) (KL) 


(1-LLT) + LLT 


(28a) 


Again, this saturation is used only when set point is more than 10 K away 
from tunnel temperature. This logic is also imposed as a hysteresis 
through LLT logic discussed below and illustrated in figure 32. 


where 

if 

ET < 

5 

, LLT = 1 


and 

if 

ET > 

10 

, LLT = 0 

(28b) 


Thus the temperature loop works on a gain scheduled pID law near the desired 
set point. When large errors are encountered, appropriate saturation and 
feed-forward compensation are used to prevent thermal shocks and energy wastage. 

Pressure loop. - Figure 31 details the schematic diagram for the control 
of tunnel total pressure. The control law used is 

(EP) = P-(SP) 

(AG) = [(pPXEP) + (IP) r (EP) + (DP) / (EP) 1 (29) 

(KG)P^ L J 
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This control law may be recognized as the control law of equation (24) . 
The pressure loop is not affected by large errors, except that when the 
tunnel pressure is low and the mach number is high, valv.e area in excess 
of 100 percent of the gas bleed valve is necessary. Hence an automatic 
law for opening two parallel analog valves as a bias has been proposed. 

This logic is 



SUM = 0 




If 

Aq > 0.9 

, SUM = 

SUM , + 
past 

0.125 

If 

Ag < 0.5 

, SUM = 

SUM 

past 

0.125 

If 

SUM < 1 

It 

1 — 1 
<d 

SUM & A2 

= 0 

If 

SUM > 1 

, A1 = 

1.0 A2 

= SUM-1 


This logic automatically opens a parallel bias valve by 12.5 percent when- 
ever the automatic control valve reaches 90 percent on control. A maximum 
of 200 percent area, is thus available at low pressures. The logic also 
automatically shuts off the bias valves. A time delay of UA is imposed 
in this logic to prevent hunting. 

Microprocessor-Based Controller 

In order to execute the complex control law on the tunnel, an electronic 
controller is necessary. The choice between an analog processor or a digital 
processor is decided by the complex nonlinear computation that needs to be 
performed. Analog controllers cannot perform such functions accurately. 

Hence, a digital microprocessor controller has been proposed. This micro- 
processor is required to continuously sense the various tunnel variables, 
compute the linear and logic components of the control law, and generate 
appropriate commands for driving the control input actuators. Following 
are the functional requirements of such a controller. 

The microprocessor is required to sense the tunnel variables, total 
temperature to a resolution of 0.1 K or better, total pressure to a resolution 
of 0.006 atm, fan speed to about 0.1 percent, fan power to about 1 percent, 
and liquid nitrogen pressure to about 1 percent. Maximum values of these 
variables are 350 K, 6 atm, 7,200 rpm, 2,000 kW and 14 atm, respectively. 
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Hence, the microprocessor should resolve the incoming signals to about 0.025 
percent resolution. This corresponds to a 12-bit AD converter having a 
6-channel capability. Typically 0- to 5-V DC signals from sensors can be 
considered. 

The bandwidth, of the tunnel temperature and pressure control systems 
is about one rad/sec. In order to implement such a control, the micro- 
processor controller must function at at least 10 cycles of computation/sec. 

In view of the complex nonlinear computations, the microprocessor arith- 
metic unit must have a floating point capability. This can be achieved only 
by floating point hardware, considering the quantity of computation to be 
performed in each cycle and the accuracy desired. 

The various I/O devices on the microprocessor controller are the AD 
converters for signal inputs, the data entry switches for the tunnel set 
point variables and the loop constants, and either direct digital or DA- 
converted analog outputs for the actuators and displays. 

For the 0.3-m TCT facility, the following I/O devices have been proposed. 
The tunnel variables, total pressure, static pressure, total temperature, fan 
speed, fan power, and liquid nitrogen pressure are digitized through a six- 
channel multiplexer and AD converter. The tunnel operator controls consist 
of dedicated decimal thumbwheel inputs for the pressure and temperature set 
points with their read buttons, and a set of data entry thumbwheels to put 
in about 60 loop constants through a key switch. The operator display consists 
of the tunnel set points, the tunnel variables, and the actuator commands, 
each for the pressure and the temperature loop. Provision for auto-hold-manual 
switching of each loop is provided. A view of the controller front panel is 
shown in figure 33. The microprocessor controller output consists of 2 10-bit 
pure binary digital valve drives for the liquid and gaseous nitrogen control 
valves, two analog outputs for liquid and gaseous nitrogen control valves 
(optionally), and two analog outputs for the bias analog valves. 

Control Flow Chart 

The detailed functional flow chart for automatic control of the 0.3-m 
TCT is provided in Appendix A. This flow chart is .'a nvnnerical mechaniza-. . 
tion of the control laws of equations (25) to (30) , shown in figures 30 
and 31 . 
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Firstly, the tunnel variables T, P, P , N, P , and the fan 
power are AD converted and scaled into appropriate engineering units. 

Then the temperature control laws of equations (25) to (28) are mechanized 
including the liquid valve gain adaptation. The p ID control law is 
invoked when in '*auto" mode. In the "manual" mode, the liquid valve is 
directly driven by the command. In the "hold" mode, the valve position is 
locked. The actuator for the liquid flow valve can be updated at a rate 
of UT cycles. 

The pressure control law of equations (29) and (30) is then mechanized 
with provision for automatic bias valve control- The pID control law is 
invoked only in the auto mode. The manual and hold modes function 
similarly to the temperature loop. 

The various statements of the flow chart in the Appendix are self- 
explanatory. These allow for necessary sensing, control laws — linear and 
logic, actuator drive, and display of 10 quantities ^ Further, provision 
has been made for fast emergency function of the pressure and the temperature 
control valves. The liquid valve fully closes on emergency, and the gas 
bleed valves open fully to avoid overpressurizing the wind tunnel. 


38 



CLOSED-LOOP CONTROL RESPONSES FOR THE 0.3-m TCT 


The control laws have been realized on a Motorola 6800 microprocessor, 
which is an 8-bit, 24 K memory processor. The control flow chart has been 
translated to an assembler language program which has been burned into 
the Read Only Memory of the processor. The high-speed arithmetic 
capability has been obtained on floating point hardware mathematical 
chips. With this unitj the microprocessor controller is capable of per- 
forming all the control functions in about 70 msec. The program is designed 
to run at a rate of 100 msec per cycle using an internal clock. 

This microprocessor controller has been successfully commissioned to 
realize pressure control to ±0,017 atm and ±0,25 K. In figure 34, a 
typical tunnel run of the 0.3-m TCT facility is shown, while under automatic 
control. Initially , a temperature set point of 100 K and; .a. pressure set 
point of 2.1 atm were sought, and microprocessor control was switched on auto 
mode. Because of large temperature- error , the LN 2 flow was limited with feed 
forward on. The tunnel pressure built up in about 30 sec. The tunnel 
temperature moved slowly to 100 K in about 30 min at a typical rate of 7 K/ 
min for the metal. Once the temperature stabilized, the pressure was 
built up to 5.5 atm. Initially a +2 K step was tried on temperature which 
was rapidly achieved. Then a number of mach number changes from 0.3 to 
0.8 were tried with a high blockage model put through many angle programs. 

The wiggles in temperature and pressure of figure 34 are attributable to 
onset of these disturbances. In every case the- temperature settled to 
within + 0.25 R,. and pressure within ±0.017 atm- of set point. During the 
tunnel run, the pressure set points were changed from 5.5 atm to 3.4 atm 
and temperature from 100 to 120 K. Throughout the run of figure 34, both 
pressure and temperature loops worked well against various disturbances. 

In figure 35, a temperature set point change and its coupling to the 
pressure loop are shown. The pressure has taken about 30 sec to stabilize 
and temperature stabilizes in about 15 sec. Since the presently used 
controller has a 10-bit ADC, a small limit cycle of about 0.25 K is occasion- 
ally evident in the temperature. Similarly, the pressure loop also has a 
small amplitude limit cycle caused by 10-bit ADC. Presently , conversion to 
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12-bit ADC is being planned and should reduce the occasional limit cycling 
to less than 0,1 K and 0.005 atm. 

In figure 36, both the pressure and temperature loops have been con- 
trolled in closed loops' to achieve - s-teppedmavement. of P and T towards 
liquid gas phase transition. This diagram shows the quality of control 
since no Mach or model attitude disturbances are present. A pressure 
control instability was encountered at very low Mach number and low 
power conditions at high temperatures. The pressure control loop functions 

on the basis of mass balance between incoming liquid nitrogen and 

« • 
gas bleed mass flow m^. Under' low-power conditions, the* m^ required 

to cancel fan heat is very low. The temperature sensitivity to is 

high at temperatures^ and hence the temperature, loop does not work fully 

in the proportional band but on the threshold of the on^roff mode. The 

pressure, control loop tracks' and hence has- a tendency for .instability 

due to discontinuous ' The problem was found, to alleviate itself as 

the power level was increased to M = 0.25. 
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CONCLUDING REMAilKS 


In this report the design of a closed-loop automatic controller 
for a closed-circuit cryogenic pressure tunnel has been detailed, and 
has been successfully used to control the 0.3-m TCT facility. The 
following observations sum up the achievements of this effort: 

(1) The mathematical control model of a cryogenic tunnel, synthe- 
sized and validated as detailed in reference 4, has been 
successfully used in the control studies, 

(2) Two single- input/single-output type control laws have been 
developed for automatic closed-loop control of the tunnel’s 
total temperature and total pressure, 

(3) The control laws devised are nonlinearly gain scheduled, 
pID control laws. These have provided satisfactory closed- 
loop control of the tunnel, firstly on the real-time cryogenic 
tunnel simulator and subsequently on the 0.3-m TCT facility, 

(4) Presently, regulation accuracies of ±0,25 K in temperature and 
±0.017 atm in pressure have been achieved on the 0.3-m TCT. 
Improvements to ±0.1 K and ±0.005 atm accuracy in regulation 
are being planned by conversion of 10-bit ADC to 16-bit ADC in 
the controller, 

C5) The control flow chart developed is a general one and can be 
used on any cryogenic pressure tunnel. 

(6) The nonlinear control logic has been found to provide a safe 
method for gross tunnel condition changes under closed-loop 
control. 

(7) Some pressure control problems were encountered at very low Mach 
number and high- temperature conditions attributable to nonlinear 
coupling from the temperature loop into the pressure loop. The 
problem corrects itself at about M = 0.25 and above. 

(8) It has been shown that the cryogenic tunnel process is basically 
well behaved and the control problem can be solved. 
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APPENDIX A 


FLOW CHART FOR CLOSED-LOOP CONTROL OF A FAN-DRIVEN 
CRYOGENIC WIND TUNNEL. 



LIST OF CONTROLS AND DISPLAYS 


CONTROLS 

!• Four digit dedicated decimal thumbwheel controls for S? 

2. Four digit dedicated decimal thumbwheel controls for ST 

3. Six digit floating/fixed point with polarity decimal thumbwheel 
controls for all other DATA 

4. Two digit dedicated decimal thumbwheel for DATA ID 

5. Push buttons (12) 

Temperature Loop: Read, Emergency, Auto, Hold, Manual 
Pressure Loop; Read, Emergency, Auto, Hold, Manual 

Other; Feed Forward, LN 2 Limit 

DISPLAYS 

1. Four 4 digit decimal displays for SP, P, ST, T 

2. Two 3 digit decimal displays for AG, ALQ 

3. One 6 digit floating /fixed point display for all DATA 

4. One 2 digit display for DATA entry ID 

5. Four 4 digit decimal displays for GFLOW, FLOW, MAX, WG 

6. LED on-off lamps for 

Two operate 
Two auto 
Two hold 
Two manual 
Two emergency 

Seven, FF on/FF off/LLT on/LLT off/MEM fault/full open/low pressure 
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Inputs Via DATA Thumbwheel Switches 


6 Channel AD Converter 

0-5VDC = Full Range Of Tunnel Variable 


Inputs Via DATA Thumbwheel 
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Mach Number 

Liquid Nitrogen Flow at Fullopen, ,kg/j 
Gain Schedule 

PID Gains 

. 

Fan Power, kV 


Temperature Error 

Fan Power Feed Forward - off 
Fan Power Feed Forward - on 
Liquid Nitrogen Flow Limit - max 
Liquid Nitrogen Flow Limit - 100% 


Bias For Feed Forward 
Liquid Nitrogen Flow Limit 


Mini - Max Limits 
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t = Time at Hold 






COMPUTE 


EP(n) = P(n) - SP 

RlP(n) = RlP(n-l) + KIP EP(n) 

IE RlP(n) 0 , RlP(n) = 0 
IE RIP(n) '^MXG , RlP(n) = MXG 
IE MXG > 1 = 1 


Pressure Loop Error 
Pressure Error Integral 

Mini-Max Limits 



Pressure Valve Drive Signal 

Integral 

Derivative 

Mini-Max Limits 
Gas Flow, kg/ sec 


HOLD 

MAMJAL 


DRIVE GAS VALVE 
AG(n) = 0 Full Close 
AG(n) = 1 Full Open 
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DIGITAL 

ANALOG 









LOGIC AND COMPUTE 

If AG(n) >0.9 SUM(n) = SUM(n-l) + 0.125 

IF AG(n) <0.5 SUM(n) = SUM (n-1) - 0.125 

IF SUM(n) I'O ; SUM(n) = 0 

IF SUM(n) ^ 1 ; A(l) = SUM(n) = 0 
A(2) = 0 

IF SUM(n) > 1 ; A(l) = 1 
A(2) = SUM(n) - 1 


DRIVE ANALOG VALVES 
A(l) = A(2) = 0 ; Full Close 


“► A(l) 

Analog Valves 

^ A(2) 


A(l) = A(2) = 1 ; Full Open 


i < UA 


i < UP 









APPENDIX B 

MACH NUMBER CONTROL LOOP 

The 0.3-m transonic cryogenic tunnel test section Mach number is 
dominantly controlled by fan speed and the tunnel gas temperature, as 
described by equation (5) of the main text: 


V — - 0.035 

T (1 - 0.3 M) P (5) 

' 

The dynamics of the tunnel test section Mach number as a function of the 
fan speed are represented by 


M 



-T 

N e 


S 

a 


k 


m 


(i 


+ t 

P 


s) 


where 


= acoustic time delay from fan to test section, sec 

t^ = plenum time constant, sec 

- 0-035 

k = 597 (1 - 0.3 M) P 
m 

The fan of the 0.3-m TCT is driven by a water-cooled, 2-pole, 3-phase 
squirrel cage motor. The motor power is derived from a variable frequency 
power source which maintains a voltage-to-frequency ratio of 35 V/Hz and 
delivers up to 2000 MW of electrical power. The fan speed is the synchronous 
speed of the motor at the given frequency less a load-dependent slip which 
is about 0.4 percent at maximiim power, and normally much less. The 0.3-m 
TCT fan speed-Mach number relationship in the subsonic range is shown in 
figure 13. The steady-state Mach number sensitivity to fan speed can be 
determined from equation C5) or figure 13. The sensitivity corresponds to 
about ±2.5 rpm for a 0.002-Mach band at 100 K and ±4.5 rpm at 300 K, The 
existing rotating speed/frequency control system is able to maintain a 
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typical speed stability of ±4 rpm, under closed-loop control. The dynamics 
of the speed control system have been determined by invoking identification 
procedures to arrive at the following description of the fan speed control 
system. 


N ^ 1 

^set 1 + 0.56 S + 0.2 $2 

This represents an average response of the fan speed control system for the 
0.3-m TCT fan load. 


Mach Number Control Loop 

The automatic Mach number control scheme for the 0.3-m TCT is shown in 
figure 37. The schematic shows the 2-pole, 35 V/Hz motor driven by a 
variable frequency generator. The fan rotation sets the cryogenic tunnel 
gas in motion to provide varying Mach numbers at the test section. The 
test section Mach number is estimated from the tunnel total pressure P and 
the test section static pressure The transducers used to sense these 

pressures have a time constant of t 2 sec. The Mach number is estimated 
as 



The Mach niomber M is compared with the set point Mach number SM. The error 
EM is gain scheduled, a design goal, and is operated on by a pi controller 
to generate the fan speed set point automatic Mach control 

mode. 

The schematic block diagram for Mach number control analysis is shown 
in figure 38. From this diagram, the open-loop gain product G(s) H(s) is 
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G(s) H(s) ^ 

^ 1 + 0.56 S + 0.2 


= T S 
a 

e 

vV* (i + v) 


(l + tas) 


= KK 


f(pM) • S + (IM)1 

S (l + 0.56 S + 0.2 $2) (l + tpS) (l + tas) 


_T S 

where e has been ignored because of the low x , which is small 

Si 

compared to t^, t 2 and other constants. The system stability is determined 
by seeking' th^ roots of polynomial 1 +*C(s) H(s) = 0 and studying their 
properties- 


1 + G(S) H(S) =0 = 0.2 ta t $5 

P 

+ [0.2 (ta -H tp) + 0.56 ta tp ] S‘* 

-H jo. 2 + ta tp - 1 - 0.56 ta tp] s3 
+ [o.56 -H ta -t- tp ] $2 
+ [l -1- (KK)(pM)] S 
-I- (IM)(KK) 

The poles of the Mach number control loop are 

s = 0; s = - ; S = - S = - 1.4 ± 1.2 j 

^ p 

The zero of the Mach number control loop is generated by the controller as 

c - (IM) 

" ■ (pM) 
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The Mach number control system has five poles and a zero resulting in four 
poles in excess of zero and hence four asymptotic trajectories. These 
asymptotes move in the direction of ^ + 45° in the S = a + ju plane. 

A sketch of the asymptotic root locus is shown in figure 38. As the 
system gain is moved from zero to a large value, two of the five poles 
move into the right half plane, indicating possible instability. By 
chposing proper gain (KK) , the Mach loop poles can be located at critical, 
damping. Since (KK) turns out to be a constant. 


KK = — = constant 

k ^ 

m 1 


KM oc k ^ 
m y 

The tunnel steady-state Mach nvimber accuracy is determined by the open-loop 
integration and hence has a zero steady-state error. The pi controller 
contributes this integration. To prevent instability associated with this 
lag, a lead-time constant of is used and is a very large number 

(typically 10 to 20 sec). 

Since the Mach loop is obviously very sensitive to large errors, the 
following saturation terms are proposed for the control law 


N 


set 


KM ((pM) • S + (IM)) 

^ (EM) 


Mini/max limits: 

If ABS (EM) >0.05 ; ABS (EM) = 0.05 

If ABS - n) > 300 rpm; ABS - n) = 300 

Since the average ^ = ^0^002^ ’ O' 05 corresponds to 150 rpm. Hence, the 
Mach loop controller limits the proportional band authority of the pi 
controller to ±300 rpm incremental authority. 
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Mach Loop Closure Implementation 

The existing speed control system for the 0.3-m TCT facility has a 
number of bands of speed singularities wherein the system cannot be operated 
for more than a short period. These singularities are caused by the rotating 
frequency converter and occur at 1000 rpm, 3500-3700 rpm, 4800-5200 rpm and 
5700 rpm. 

A simpler Mach loop' closure^ based on the. pi .control law and simple, 
error-dependent linear/logic control laws are not intelligent enough to 
avoid these bands for a certain set of specific Mach numbers. Hence, it is 
considered inappropriate to close the Mach loop unless corrections are 
made for the speed singularities. 

Presently used manual speed set point control of the tunnel provides 
adequate manual control. 


Concluding Remarks 

In this appendix, the design of a pi controller for the Mach number 
control loop has been discussed. The gain scheduling and pi settings 
necessary for stable Mach number control have been discussed. The bandwidth 
of the Mach loop is about 1.2 rad/sec and has a small perturbation settling 
time of 15 sec. 
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MACH LOOP CONTROL FLOW CHART 


CLOCK CONTROL^ 
START 


AD 

DEL 

S3 

Convert and scale into engineering units 
Sampling time, Sec 

n 

= 

Time index, integer 

P (n) 

= 

Total pressure. Atm 

PS(n) 

= 

Test section static. Atm 

T(n) 


Total temperature 

N(n) 

= 

Fan speed, rpm 


M = 


ESTIMATE 

Measured Mach number 



Analog Inputs 


0-5VDC 

0-5VDC 

0-5VDC 

0-5VDC 


0-6.8 atm 
0-6.8 atm 
75-325 K 
0-6400 rpa 
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Mach Loop Error 



Panel Controls : 

4 digit displays each for M(n) , NS(n), SM 

Push buttons for Auto, Hold, Manual, Emergency, Read 

Thumbwheel switches 4 decimal digit for Set Point 
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Figure 1. 0.3-m TCT control configuration. 
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Figure 2. Liquid nitrogen charge and gaseous nitrogen bleed schematic. 
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5. Tunnel average gas time constant. 
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Assumptions 

Perfect gas behavior 

Uniform tunnel temperature 

Work, potential & kinetic energy ignored 


first law of thermodynamics 




W C -r 
g V dt 


u(m. - m ) 

L g 


Where 

h^^ Liquid nitrogen enthalpy J/kg 
u Specific internal energy 
hg Gas enthalpy J/kg 

a Heat flow from tunnel metal wall J/sec 


Qp Heat flow from Fan operation j/sec 
m Gaseous nitrogen bleed kg/sec 

g 

Liquid nitrogen mass flow kg/sec 


Figure 6. Thermodynamic model of a cryogenic tunnel. 
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Figure 7. Multivariable model of a cryogenic tunnel. 




Temperature K 


Figure 8. Steady-state temperature gain for LN 2 input. 
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Temperature K 


Figure 10. Steady-state temperature gain for GN 2 bleed. 
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ACTUATORS AND SENSOR MODELS 


LN^ VALVE - = k^A^^ = 3. 47 - P A^_ 


GN, VALVE - A = 21. 8 -p A P > 1. 5 ' 

2 g g VI g Vt 9 


= 21. 8 

TEMPERATURE SENSOR 
PRESSURE SENSOR 
FAN DYNAMICS 


V"T 

1 


2 -• 


U 

L P J 


1.7 


I 


I 


1 + t^s 


N 


A KP<1.5 

y g - - 

WHERE t = - TO 3 secs 
1 2 

WHERE L = isecs 
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5 
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N 
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Figure 14. Actuator and sensor dynamic models for 0.3-m TCT. 
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GAS VALVE 



Figure 15. Cryogenic tunnel temperature control schematic. 
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Figure 16. Asymptotic root loci for temperature loop. 
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Figure 17. Root loci for temperature loop — Mach number effects. 
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Figure 18. Root loci for temperature loop — temperature effects. 


77 



LIQUID TEMPERATURE 
VALVE LOOP 


PRESSURE 
SET POINT 





j. 


Figure 19. Cryogenic tunnel pressure control schematic. 
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Figure 20. Asymptotic root locus for pressure loop. 
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Figure 21. Noninteractive multivariable controller. 




Figure 22. View of the hybrid computer-based cryogenic t-nnel simulator. 
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Figure 24. Tunnel temperature dynamics — analog. 
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Figure 26. Tunnel valve control scheme — analog. 
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Figure 27. Tunnel pi controllers — analog. 
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PRESSURE SEHLING TEMP. SETTLING 

TIME ^ 20 SEC TIME 20 SEC, 



Figure 28.. Simulator closed-loop control responses (270 K) . 
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PRESSURE SETTLING TEMP. SETTLING 
TIME** 20 SEC TIME* 20 SEC 



Figure 29. Simulator closed-loop control responses (90 K) . 
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Figure 30. Automatic temperature control scheme for a cryogenic tunnel. 
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i 


Control law 

EP • P- S? 

r 1 

" p pP'EP + Ipfpp OP* i cp 
kg p^ 7 L ot 


logic 

SUM -0 ^^^past ’ seconds ago 

IF AG>0,9 : SUM • SUM . +0,125 
past 

IF AG<a5 : SUM-SUMpgjj -0.125 
IF SUM<1 ; A I • SUM A 2 -0 

IF SUM>1 : A 1 • UO A 2 • SUM -I 


Figure 31. Automatic pressure control scheme for a cryogenic tunnel 






Fan power feed forward 


On 



Temperature error, K 


Liquid nitrogen flow limit 



Temperature error , K 


Figure 32. Cryoger 'c tunnel temperature control loop logic functions. 
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Figure 33. View of the microprocessor-based controller for 0.3-m TCT. 
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Figure 34. Typical 0.3-m TCT test with temperature and pressure automatic 
control. 
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Figure 35. Typical temperature response and pressure coupling under closed- 
loop control. 
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Figure 36. Typical 0.3-m TCT test with temperature and pressure automatic 
control . 
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Figure 37. Automatic Mach number control scheme for a cryogenic wind 
tunnel . 
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Figure 38. Cryogenic wind tunnel Mach number control analysis. 
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